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LONG-TERM  GOALS 

Modeling  and  predicting  ocean  optical  properties  for  coastal  waters  requires  linking  optical  properties 
with  the  physical,  chemical,  and  biological  processes  in  the  upper  ocean.  Our  long-term  goal  is  to 
incorporate  optical  processes  into  coupled  physical-biological  models  for  coastal  waters,  develop  and 
improve  integrated  ocean  forecasting  systems,  including  prediction  of  ocean  optical  properties. 

OBJECTIVES 

1)  To  improve  performance  of  the  coupled  physical-biological  model,  which  is  based  on  the  Navy 
Coastal  Ocean  Model  (NCOM)  and  Regional  Ocean  Model  System  (ROMS)  for  the  California  Current 
System  (CCS); 

2)  To  incorporate  optical  variables  into  the  improved  coupled  3D  physical-biological  model  for  the 
CCS; 

3)  To  use  these  variables  to  drive  a  radiative  transfer  model  (EcoLight)  that  will  simulate  and  predict 
the  subsurface  light  field  as  well  as  the  ocean’s  color; 

4)  To  conduct  NCOM-biological-optical  model  simulations  for  the  Monterey  Bay  for  the  period  of  the 
intensive  field  experiments,  such  as  BioSpace  in  Monterey  Bay  2008  (June  2008),  and  use  the  bio- 
optical  measurements  to  constrain  and  improve  the  model  performance. 

APPROACH 

To  achieve  the  first  objective,  we  have  conducted  a  series  of  3D  physial-biological  model  simulations 
for  the  California  Current  System.  We  have  incorporated  the  Carbon,  Silicate,  and  Nitrogen  Ecosystem 
(CoSiNE)  model  into  both  the  Navy  Coastal  Ocean  Model  (NCOM)  and  Regional  Ocean  Model 
System  (ROMS).  The  CoSiNE  model  (Chai  et  ah,  2002;  2003,  and  2007)  was  developed  originally  for 
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the  equatorial  and  Pacific  Ocean.  During  the  past  three  years,  we  have  implemented  the  CoSiNE  model 
and  improved  the  perfonnance  for  the  CCS,  and  investigated  the  seasonal  and  interannual  variation  of 
biological  productivity  in  response  to  both  local  and  remote  climate  forcing. 

For  the  optical  modeling  work,  our  approach  is  to  incorporate  spectrally-resolved  inherent  optical 
properties  into  the  existing  NCOM-CoSiNE  and  ROMS-CoSiNE  model  which  haven’t  included  the 
optical  processes  and  in-water  radiative  transfer  model  for  the  California  Current  System  (see  Figure 
1).  We  have  successfully  incorporated  an  optical  module  along  with  a  commercially  available  radiative 
transfer  model  (EcoLight)  into  a  one-dimensional  ecosystem  model  which  is  flexible  enough  that  can 
be  applied  to  three-dimensional  circulation  model  (Fujii  et  ah,  2007).  This  effort  illustrates  the  value  of 
adding  optical  processes  into  physical-biological  models.  Based  on  this  experience,  we  developed  an 
optical  model  that  is  fully  compatible  with  NCOM-CoSiNE  model  for  the  California  Current  System. 
The  optical  model  explicitly  represents  spectrally-resolved  inherent  optical  properties  (IOPs)  including 
absorption,  scattering,  and  attenuation.  The  optical  model  requires  input  from  NCOM-CoSiNE  model, 
such  as  phytoplankton  associated  chlorophyll  and  non-algal  particles  (NAP),  and  it  provides  outputs 
such  as  phytoplankton  absorption  coefficient,  detritus  absorption  coefficient,  particulate  backscattering 
coefficient.  These  optical  products  are  all  spectrally-resolved  and  can  be  compared  with  satellite  and  in 
situ  data  directly. 


Figure  1:  Overall  approach  for  incorporating  optical  properties  into  NCOM-CoSiNE  model. 
The  effect  of  CDOM  in  the  optical  model  is  a  constant  at  present,  but  will  be 
modeled  dynamically  in  the  next  phase  of  the  project. 


In  the  optical  model,  absorption  coefficient  is  determined  from  the  sum  of  the  absorption  coefficients 
of  seawater,  phytoplankton  (based  on  the  chlorophyll  content  of  pico-plankton  and  diatoms, 
respectively),  NAP,  and  CDOM.  Seawater  absorption  coefficient  is  modeled  based  on  previous 
observations  (Pope  and  Fry,  1997;  Morel,  1974).  Phytoplankton  absorption  is  calculated  as  the  sum  of 
pico-plankton  and  diatoms  with  different  chlorophyll-specific  absorption  coefficients  due  to  the 
different  sizes  of  them.  Pico-plankton  has  a  higher  chlorophyll-specific  absorption  coefficient  and  a 
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steeper  absorption  spectra  than  diatoms.  In  addition,  the  absorption  spectra  of  a  given  phytoplankton 
functional  group  changes  with  intercellular  chlorophyll  concentration  (the  package  effect).  Therefore, 
the  chlorophyll-specific  absorption  coefficients  by  pico-plankton  and  diatoms  are  modeled  separately: 

a^(A,z)  =  a*JA,z)xchl\(z)  +  a*l/l2(A,z)xchl2(z)  (1) 

where  chll  and  chl2  are  chlorophyll  concentration  for  pico-plankton  and  diatoms  converted  from 
modeled  biomass.  A  photo-acclimation  model  is  included  to  include  the  light  adaptation  effects  on 
chlorophyll/carbon  ratio. 

Color  dissolved  organic  matter  (CDOM)  also  plays  important  role  in  absorption  spectra.  We  fix  its 
absorption  coefficient  at  440nm  vertically  to  be  0.016  (m'1),  which  is  within  observed  range  (0.012- 
0.019  (nT1),  Simeon  et  ah,  2003).  We  do  so,  as  the  first  step,  because  the  current  version  of  NCOM- 
CoSiNE  model  has  not  included  active  CDOM  dynamics  and  ratios  of  CDOM  to  dissolved  organic 
matter  (DOM)  or  dissolved  organic  carbon  (DOC),  and  these  processes  vary  greatly  in  the  CCS.  In  the 
optical  model,  CDOM  absorption  spectra  is  assumed  constant  vertically  and  with  a  spectral 
dependence: 

aCDOM (2)  =  0.016  x  exp{-0.0145  x  (A  -  440)}  (2) 

NAP  absorption  follows  a  similar  manner  to  CDOM,  and  is  expressed  as: 

aNAP(2,z )  =  a*NAP{ 440)  x  NAP(z)  x  exp {-0.01 1  x  (T  -  440)}  (3) 

where  aNAp+(440)  is  the  carbon-specific  absorption  coefficient  by  NAP  at  440nm  which  is  assumed  to 
2  1 

be  0. 1  (m  gC"  )  based  on  observations  (Babin  et  al.,  2003). 

In  the  optical  model,  assuming  no  contribution  of  CDOM  to  backscattering,  backscattering  coefficient 
is  determined  as  the  sum  of  seawater,  algae  and  co-varying  particles.  Algae  and  co-varying  particles 
are  expressed  as  a  function  of  particulate  organic  carbon  concentration  (POC)  of  small  and  large 
particles,  which  consist  of  pico-plankton,  diatoms,  and  associated  NAP,  respectively.  Backscattering 
by  small  and  large  POC  (POC1  and  POC2)  are  formulated  as  follows  (Stramski  et  al.,  1999): 


h  (T,z)  =  (P°a(z))i-277  x(  2  )-°-5 
bpl  476935.8  510 

(4) 

b  (A  z)  =  (POC2^z\oAs9 
bpA  ,  )  (  17069.0  j 

(5) 

The  total  backscattering  coefficient  is  calculated 

bh(A,z)  =  0.5  bsw(A)  +  bbpl(A,z)  +  bhp2(A,z)  +  bbbg  (6) 
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where  bsw(A)  is  seawater  scattering  coefficient  with  the  correction  for  salts  (Boss  and  Pegau,  2001). 
bbbg  is  the  background  backscattering  coefficient  the  implicitly  reflects  contributions  by  non- 
phytoplankton-covarying  bacteria  and  other  particles  (Behrenfeld  et  ah,  2005). 

WORK  COMPLETED 

This  report  summarizes  modeling  activities  between  1  January  2008  and  30  September  2008. 

We  have  completed  the  following  modeling  tasks: 

1 .  Evaluated  the  performance  of  3D  physical-ecosystem  model  for  the  California  Current  System 
(CCS)  with  the  observations  from  the  Californian  Cooperative  Ocean  and  Fisheries 
Investigation  (CalCOFI)  and  Monterey  Bay  Aquarium  Research  Institute  (MBARI); 

2.  Finked  the  bio-optical  module,  based  on  ID  configuration  (Fujii  et  al.,  2007),  with  the  Carbon, 
Silicate,  Nitrogen,  and  Ecosystem  (CoSiNE)  model  for  the  CCS; 

3.  Incorporated  bio-optical  module  and  CoSiNE  into  the  Navy  Coastal  Ocean  Model  (NCOM)  for 
the  California  Current  System  (CCS),  and  the  Regional  Ocean  Model  System  (ROMS)  for  the 
Pacific  Ocean; 

4.  Conduct  a  series  of  model  simulations  to  test  the  robustness  of  the  NCOM-CoSiNE-Optics  and 
ROMS-CoSiNE-Optics  for  various  regions; 

5.  Started  to  evaluate  model  perfonnance  with  remote  sensing  derived  IOPs  for  the  CCS. 

RESULTS 

We  have  conducted  a  series  of  the  ROMS-CoSiNE  model,  without  the  optical  component,  for  the 
California  Current  System  for  the  period  of  1990-2005.  For  doing  so,  we  can  evaluate  the  ROMS- 
CoSiNE  model  performance  with  the  long  term  time  series  observations  from  the  Californian 
Cooperative  Ocean  and  Fisheries  Investigation  (CalCOFI)  and  Monterey  Bay  Aquarium  Research 
Institute  (MBARI).  First,  compare  the  modeled  physical  and  biological  variables  on  the  annual  mean 
basis,  then  evaluate  the  seasonal  cycle  of  nutrient  dynamics  and  biological  productivity.  The  CalCOFI 
and  MBARI  observations  are  grouped  into  six  dynamically  different  regions  according  to  near-shore 
upwelling  areas  vs.  off-shore  less  productive  regions,  Figure  2. 
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Figure  2:  Long  term  CalCOFI  and  MBARI  observations  for  the  central  California  regions.  Up  to 
166,612  profiles  of  measurements  from  CTD  casts  (temperature,  salinity,  nutrients,  chlorophyll, 
oxygen)  were  available  as  well  as  1519  records  of  primary  productivity  measurements  and  886 
records  of  taxonomic  composition  for  phytoplankton  in  the  Monterey  Bay. 


The  observed  data  and  the  modeled  simulated  results  are  compared  for  each  sub-regions,  and  the 
statistical  analysis  has  been  performed.  Overall,  the  model  behavior,  in  terms  of  chlorophyll  and  nitrate 
concentrations,  is  quantitatively  consistent  with  observations  (high  r2)  especially  for  the  off-shore 
domains.  Figure  3  shows  the  annual  mean  modeled  surface  chlorophyll  and  nitrate  for  the  entire 
domain.  The  ROMS-CoSiNE  model  captures  the  on-  and  off-shore  contrasts,  is  able  to  reproduce  the 
spatial  extension  of  the  upwelling  system. 
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Figure  3:  Annual  mean  modeled  surface  chlorophyll  and  nitrate  comparisons  with  the  long-term 
observations  (CalCOFI  and  MBARI),  Pawlowski  et  al.  (2008). 

Beside  the  incorporation  of  the  CoSiNE  model  into  the  ROMS,  we  have  been  collaborating  with  the 
Dynamics  of  Coupled  Processes  Section  (led  by  Dr.  Igor  Shulman)  at  NRL  to  implement  our  CoSINE 
model  into  the  Navy  Coastal  Ocean  Model  (NCOM)  for  the  California  Current  System  (CCS).  The 
coupled  physical  and  biological  model  assimilates  remote  sensing  information,  and  produces  near-real 
time  simulations  on  regular  basis,  see  Figure  4.  The  physical-biological  model  results  can  be  viewed 
and  obtained  directly  from  the  NRL  website: 
http://www7320.nrlssc.navy.mil/ccsnrt/ncomccs-9Km_main.html 
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MODEL  CHLOROPHYLL  AND  CURRENTS 
25-SEP-2008  to  29-SEP-2008 


— L3D  -128  -126  -124  -122  -120  -US  -116 


Figure  4:  The  NCOM-CoSiNE  modeled  surface  currents 
and  chlorophyll  for  the  CCS,  25-29  September  2008. 


We  use  this  NCOM-CoSiNE  model  as  our  baseline  model  configuration,  and  incorporated  the  optical 
module  into  this  near  real-time  modeling  system  for  the  CCS.  We  used  the  NCOM-CoSiNE  model 
output  (chlorophyll,  NAP,  and  POC)  for  the  summer  of  2008  (from  May  to  August)  to  drive  an  off-line 
optical  output  (equations  1  to  6  in  this  report,  and  Figure  1).  We  evaluate  the  optical  properties  with  the 
remote  sensing  derived  inherent  optical  properties  (IOPs)  based  on  the  quasi-analytical  algorithm 
(QAA)  (Lee  et  ah,  2002),  Figure  5.  Spatial  averaged  (for  the  CCS)  total  absorption  a(440)  is  0.0489 
(m'1)  from  the  model  results,  and  0.0939  (m"1)  for  the  MODIS  data  (based  on  QAA);  the  phytoplankton 
absorption  a(p(440)  is  0.0264  (m'1)  for  the  model  and  0.0268  (m'1)  for  the  MODIS  data;  the  total 
backscattering  coefficient  bb(550)  is  0.0038  (m'1)  for  the  model  and  0.0045  (m'1)  for  the  MODIS  data. 
Overall,  the  optical  model  produced  the  IOPs  tend  to  agree  well  with  the  MODIS  derived  results  in 
terms  of  the  magnitude  and  spatial  patterns.  (Note  that  ocean  color  inversion  parameters  have  an 
inherent  uncertainty  of  about  50%,  IOCCG  report  5  edited  by  Lee).  We  are  analyzing  more  model 
results,  and  improving  the  optical  model  performance. 
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Figure  5:  Monthly  (June  2008)  averaged  comparison  between  the  modeled  and  QAA  derived  total 
absorption  coefficient  at  440  nm,  a(440)  (the  top  panel);  the  phytoplankton  absorption  coefficient 
at/440)  (the  middle  panel);  and  total  backscattering  coefficient  at  550  nm  b/,(550)  (the  lower  panel). 

The  left  column  is  the  NCOM-CoSINE-optics  model  outputs,  the  central  column  is  MODIS  QAA 
products,  and  the  right  column  is  the  statistical  comparison.  The  color  bar  is  for  the  left  and  central 
columns,  has  the  unit  of  in1.  The  colors  in  the  statistical  analysis  (the  right  column)  denote  the 

density  distribution. 


Our  current  NCOM-CoSiNE  model  configuration  hasn’t  included  the  dynamical  variation  of  CDOM  in 
light  absorption,  instead  used  a  constant  absorption  coefficient  (the  equation  2).  During  the  next  year 
of  this  project,  we  plan  to  add  dissolved  organic  matter  (DOM)  dynamics  to  the  CoSiNE  model  (note: 
we  have  finished  that  in  one-dimensional  CoSiNE  model  configuration,  which  is  flexible  enough  to  be 
applied  to  the  3D  NCOM  and/or  ROMS).  We  will  investigate  the  relation  of  CDOM  and  DOM  in 
California  Current  System  focusing  on  Monterey  Bay  using  observational  data  to  assess  how  the 
CDOM-DOM  relation  changes  in  space  and  time.  We  plan  to  use  the  MODIS  derived  CDOM 
absorption  values  to  guide  and  constrain  the  CDOM  performance  in  our  model,  Figure  6. 
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Figure  6:  Monthly  (June  2008)  averaged  MODIS  QAA  derived  surface  CDOM  and  detritus 
absorption  coefficient  at  412nm  ( acdm(412 ))  for  the  California  Current  System.  The  CDOM 
absorption  is  much  larger  than  the  detritus  component  in  this  estimate. 


IMPACT/APPLICATIONS 

Incorporating  ocean  optical  processes  into  coupled  physical-biological  models  will  enable  us  to 
simulate  and  forecast  optical  properties  in  coastal  waters.  With  demonstration  of  some  initial  successes 
of  developing  physical-biological-optical  modeling  and  data  assimilation  capability  for  the  California 
Current  System,  we  should  be  able  to  start  the  development  of  an  end-to-end  ocean  forecasting  system. 
Such  modeling  system  would  be  a  powerful  tool  to  design  the  adaptive  sampling  strategy  and  would  be 
an  essential  component  of  future  field  experiments  both  in  and  outside  Monterey  Bay. 

TRANSITIONS 

We  have  been  working  with  Drs.  John  Kindle,  Igor  Shulman,  Bred  Penta,  and  Sergio  Derada  at  the 
Naval  Research  Laboratory  (NRL).  The  ecosystem  model  code  has  been  transferred  to  the  Dynamics 
of  Coupled  Processes  group  (led  by  Dr.  Igor  Shulman)  at  the  NRL.  Once  the  optical  module  is  fully 
tested  with  the  CoSiNE  for  the  CCS,  then  the  optical  module  along  with  the  updated  versio  fo  CoSiNE 
model  will  be  transferred  to  NRL.  Fei  Chai  has  scheduled  to  travel  to  the  NRL  in  early  September  of 
2008  to  work  on  some  issues  regarding  to  the  ecosystem  component  and  the  optical  module  in  the 
NCOM,  but  we  had  to  cancel  the  meeting  due  the  impact  of  the  Hurricane  Gustav.  Fei  Chai  will  visit 
NRL  at  the  end  of  2008.  We  have  regular  communications  between  the  NRL  and  the  University  of 
Maine  on  this  collaboration. 

RELATED  PROJECTS 

This  project  has  strong  collaboration  with  other  ONR  supported  projects.  Besides  working  closely  with 
the  modeling  group  at  the  NRL  and  their  BioSpace  project,  we  are  collaborating  with  Dr.  Curtis 
Mobley  of  Sequoia  Scientific  Inc.  on  improving  the  link  between  the  radiative  transfer  model 
(EcoLight)  within  the  NCOM-CoSiNE-Optics.  Recently,  Dr.  Mobley  has  obtained  a  grant  from  ONR 
to  improve  the  efficiency  of  the  EcoLight  calculation  in  the  3D  framework.  We  are  also  collaborating 
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with  scientists  at  the  Monterey  Bay  Aquarium  Research  Institute  (MBARI)  to  use  the  observational 
data  for  the  region.  Dr.  Yi  Chao  at  JPL  has  been  collaborating  with  us  about  implementing  the  CoSiNE 
into  the  ROMS  for  the  Pacific  Ocean  and  the  CCS. 
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